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Abstract 

Shear relaxation and dynamic density fluctuations are studied within a Rouse 
model, generalized to include the effects of permanent random crosslinks. We 
derive an exact correspondence between the static shear viscosity and the 
resistance of a random resistor network. This relation allows us to compute 
the static shear viscosity exactly for uncorrelated crosslinks. For more general 
percolation models, which are amenable to a scaling description, it yields the 
scaling relation k = (p — j3 for the critical exponent of the shear viscosity. Here 
(3 is the thermal exponent for the gel fraction and (j) is the crossover exponent 
of the resistor network. The results on the shear viscosity are also used in 
deriving upper and lower bounds on the incoherent scattering function in the 
long-time limit, thereby corroborating previous results. 
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I. INTRODUCTION 



The most prominent critical phenomena at the gelation or vulcanization transition are 
of dynamic origin, see e.g. Ref. |l] for a recent review. Relaxation times in the sol phase 
diverge as the gelation transition is approached,! giving rise to critical behaviour of the 
transport coefficients: The effective diffusion constant goes to zero and the static shear 
viscosity diverges,ili indicating structural arrest of a macroscopic fraction of the monomers; 
see Figure |1] below for a graphical illustration of the situation. There are strong precursors of 
the gelation transition in the fluid-like sol phase: The decay of correlations is not exponential, 
but follows a Kohlrausch law.li This indicates a spectrum of relaxation times extending to 
arbitrary large values. 

All these phenomena are reminiscent of the glass transition so that the gelation transition 
can be regarded as a paragon for the latter. There are however important differences which 
are best explained in the context of vulcanization. In that case the distinction between 
thermal and quenched degrees of freedom is simple and unambiguous: Chemical crosslinks 
are quenched and the monomers equilibrate in a fixed crosslink configuration. In contrast to 
this scenario, the quench in the structural glass by lowering the temperature does not allow 
for such a simple classification. Even though it is generally believed that while decreasing 
the temperature at a finite rate uj some degrees of freedom are quenched, namely those 
whose relaxation times r are larger than the inverse quench rate, tuj > 1. Nevertheless 
the identification of these quenched degrees of freedom is neither simple nor unique, and 
some properties of structural glasses do in fact depend on the rate of quench.0 The so-called 
physical gelation seems a better candidate to model glassy behaviour. In that case the tem- 
perature can be easily adjusted so that the binding energy for two monomers participating 
in a crosslink is comparable to the thermal energy. Crosslinks form and break up so that 
both monomer positions and crosslinks equilibrate. By lowering the temperature crosslinks 
are quenched and give rise to structural arrest at sufficiently low temperatures. 

In this paper we will concentrate on vulcanization and chemical gelation, both of which 
are commonly interpreted as a percolation transition.i The transition from a fiuid phase to 
an amorphous solid phase happens when a macroscopic cluster of crosslinked polymers has 
been formed. However, percolation theory can only account for the geometric connectivity of 
the macromolecules and neither thermal fiuctuations nor dynamic phenomena are comprised 
in the percolation picture. Instead, one should start from a dynamic model, as is done here. 
The simplest model. Rouse dynamics, ignores all interactions except for those which ensure 
the connectivity of the clusters. De Genne^ was the first to estimate the static shear 
viscosity near the gelation transition. Relating the viscosity rjin) of a cluster of mass n to 
the longest relaxation time of the cluster, he argues that rj{n) ~ scales like the squared 
linear dimension of the cluster. The radius of gyration is determined by the mass of the 
cluster according to i? ~ n^^^f , where df = d — P/u is the Hausdorff-Besicovitch dimension 
of the fractal and related to the exponent for the correlation length u and the gel fraction f3. 
Given a scaling ansatzS for the number of n-clusters per polymer, r„ ~ n~'^ exp{—n/n*} with 
^* r-u one can average over all cluster sizes to obtain r] = '^n'^{n)nTn ~ \p—Pc\~^ 

with k = 2v — (3. This scaling relation was rederived in many other ways. It is incorrect 
because it ignores the internal structure of percolation clusters which is not only determined 
by the exponents /3 and v. The latter rule the behaviour of clusters on large spatial scales. 
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FIG. 1. Sketch of some critical quantities for gelation as a function of the crosslink concentra- 
tion c for a mean-field distribution of crosslinks. The gelation transition occurs at Ccrit = 1/2. The 
plot displays the averaged static shear viscosity (solid line) according to Eq. ( ^.261) , the effective 
diffusion constant (dashed line) according to Eq. ( |4.17| ) and the gel fraction (dotted line) according 
to Ref. g 

whereas the intrinsic self-similar connectivity is characterized by the spectral dimension d^. 
In contrast to a linear chain with (i^ = 1, a percolation cluster contains a hierarchy of 
branches and possibly loops, giving rise to a value oi dg > 1. Apparently, Catell was the 
first to notice that the relaxation spectrum acquires a Lifshits tail which modifies the above 
estimate for the viscosity //(n) of a single percolation cluster. 

In a different approach the viscosity at the sol-gel transition has been related to the con- 
ductivity of a random resistor networkcj which is made up of a fraction p of superconducting 
and a fraction [1 — p) of normal conducting bonds. Approaching the percolation transition 
from below, the conductance a ~ \p — Pc\~'^ is expected to diverge with an exponent s, 
which was predicted to be the same as k. Before it was realized that percolation clusters are 
multifractals, the two scaling arguments were considered to give identical results, namely 
k = s = 2v — (3. If multifractality is treated properly, the two scaling arguments give dif- 
ferent results. This is most easily seen in two dimensions, where dualityliil implies s = 4>. 
Here, is the so-called crossover exponent, which cannot be expressed through j3 and u. It 
was computed by Harris and LubenskyEl in the context of random resistor networks where 
a fraction 1 — p of normal conducting bonds has been removed. The crossover exponent (f) 
characterizes the average resistance il{r) ~ r*^/^ between two connected nodes which are a 
distance r apart. 

In this paper we concentrate on the fluid-like sol phase of a gelling polymer melt. We 
show that for Rouse dynamics the exponent of the static shear viscosity k is given exactly 
by 

k = (j)-(3 (1.1) 

in disagreement with both of the above scaling arguments. The derivation of the exponent 
k requires two steps. First, the viscosity is expressed in terms of the random connectivity 
matrix, which characterizes the connections in the polymer network. In a second step we 
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prove that bonds between monomers can be identified with electrical resistors and use results 
from the theory of random graphs to relate matrix elements of the connectivity matrix to 
the resistance between two nodes i and j of the corresponding random resistor network. In 
fact, we consider the derivation of such a correspondence as one of the main results of this 
paper. 

For a mean-field distribution of crosslinks we compute the static shear viscosity ex- 
actly and find a logarithmic divergence, in agreement with the general scaling relation 
( |1 . 1| ) for mean- field exponents or space dimensionality d > 6. For a crosslink distri- 
bution which corresponds to finite-dimensional percolation, we know the e-expansion of 
k = e/6 + lle^/1764+(9(e^) or use high-precision simulations in = 3 to obtain k\d=3 ~ 0.71. 

The established correspondence between the random polymer network and the random 
resistor network turns out to be useful also for the computation of the intermediate incoher- 
ent scattering function St{q). The long-time asymptotics of St{q) was discussed in Ref. ||for 
a mean-field distribution of crosslinks. In this paper, we derive upper and lower bounds on 
St{q) which put our previous results on a firmer basis. We furthermore use a scaling analysis 
to extract the long-time asymptotics of St{q) for a crosslink distribution which corresponds 
to finite-dimensional percolation. 

The paper is organized as follows: In the following Section ^ we introduce the basic 
dynamic model of a monodisperse sol of phantom monomer chains, which are permanently 
crosslinked by Hookean springs chosen at random.lli We present the formal solution of the 
dynamic model and introduce the two observables, which we want to discuss: the stress 
relaxation function and the incoherent scattering function. In Section |T| we discuss the 
static shear viscosity. First, the relation between our dynamic model and random resis- 
tor networks is established. Second, we relate the critical exponent of the viscosity to the 
exponents arising in a scaling description of the crosslink distribution, as exists for finite- 
dimensional percolation. Third, we show how to compute the static shear viscosity exactly 
for the case where the crosslinks are distributed according to mean-field percolation. This 
result is then rederived using replicas in the last subsection of Section |T|. The intermediate 
incoherent scattering function is discussed in Section IV. We first derive exact upper and 
lower bounds, which are then used to extract the long-time asymptotics for a mean-field 
distribution of crosslinks and corroborate previous results.! Some details of this calculation 
are delegated to the Appendix. In the last part of Section |rV| we compute critical exponents 
associated with the long-time asymptotics of the scattering function for a crosslink distri- 
bution which admits a scaling description — as is the case for finite-dimensional percolation. 
We conclude with some remarks and perspectives for future directions in Section ^ A brief 
account of our results was given in Ref. p. 



II. DYNAMIC MODEL AND ITS SOLUTION 
A. Dynamic Model 

We consider a system of N identical, mono-disperse (macro-) molecular units, each con- 
sisting of L monomers. These units may be monomers, dimers, linear chains, or arbitrary 
branched structures, which may also contain loops. In the following these molecular units are 
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called polymers, regardless of their topology or their degree of polymerization L. Monomer 
s on polymer i is characterized by its time-dependent position vector Ht{i, s) {i = 1, . . . , N 
and s = 1, . . . , L) in (i-dimensional space. We suppose that the polymers are subjected to a 
space- and time-dependent external velocity field (r). Here, Greek indices indicate Carte- 
sian co-ordinates a = x,y, z, . . ., and we will always consider a flow field in the r^-direction, 
increasing linearly with r^, i.e. 

<(r) := 5a,./^^r^ (2.1) 

with a time-dependent shear rate Hf As usual, 6a,i3 denotes the Kronecker symbol, i.e. 
^a,i3 = 1 for a = j3 and zero otherwise. 

We employ the simplest, purely relaxational dynamics, see e.g. Chap. 4 in Ref. [15| or 
Chap. 15.1 in Ref. [T| 

1 (^J-f 

dtRtih s) = -- + <(R,(^, s)) + er(z, s). (2.2) 

In the absence of the thermal noise ^ the monomers relax to the state dH/dH. = such 
that their velocity is equal to the externally imposed velocity field. The noise ^ has zero 

mean and covariance (i, s) s') = 2(~^ Sa,i3 Si^i' 6s,s' S{t — t'), where S{t) is the Dirac 5- 
function. Here, the overbar indicates the average over the realizations of the Gaussian noise 
^. The relaxation constant is denoted by (, and we use energy units such that ksT = 1. 



In (2^) the Hamiltonian H := Hw + U consists of two terms. The first one guarantees the 
connectivity of each polymer, whose internal structure is characterized by a connectivity 
matrix rpoiy(s; s'), 

^ N L 

^^■=^Y.Y1 rpoiy(s; s') R(z, s) ■ R(z, s') . (2.3) 

i=l s,s'=l 

We have chosen a harmonic potential to constrain the relative distance between two 
monomers on the same polymer with the typical distance given by the persistence length 
/ > 0. As an example for Tpoiy we mention the special case of linear chains. Their connec- 
tivity matrix is explicitly given by 

L-l 

rpoly(s; S') = '^{Ss,r - 5s,r+l) (^s'.r " 5s',r+l) • (2.4) 
r=l 

The second part of the Hamiltonian models M permanently formed crosslinks between 
randomly chosen pairs of monomers. The configuration of crosslinks will be denoted by 
Q = {ie, Se^i'e, s'^}^Li, i.e. the list of the M crosslinked pairs of monomers. The relative 
distances between any two monomers, participating in a crosslink, are constrained by a 
harmonic potential 

, M 

U■.= —,J2(K{^e,s^)-K{^'^,s'^))^ (2.5) 

e=l 
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whose strength is controlled by the parameter a > 0. For a ^ hard crosslinks can 
be recovered.lli Since the Hamiltonian is quadratic in the monomer positions, it can be 
expressed in terms of a NL x NL connectivity matrix F according to 

, N L 

^ =^ ^ E E s- 1\ s') R(z, s) ■ s') . (2.6) 

i,i'=l s,s'=l 

The connectivity matrix F has a deterministic part Fpoiy, which reflects the internal structure 
of the polymers, and a random part representing the crosslinks 

9 M 

F(i, s] i', s') = — 5j,i'Fpoiy(s, s') + ^(5i,j^(5s,se - 5i,i'^^s,s'^){5i',i^5s',s^ - Si'^i'^6s',s'J ■ (2.7) 

In order to determine the model completely, it only remains to specify the probability 
distribution of the crosslink configurations. We shall discuss two different types of probability 
distributions P{Q): 

i) crosslinks are chosen independently with eaual probability for every pair of monomers, 
corresponding to a mean-field distribution,0'li3 

ii) a distribution of crosslinks, which generates clusters amenable to the scaling description 
of finite-dimensional percolation. 

The precise characterization of these distributions is given below. 



B. Observables 

We shall focus on two observables. First, we aim at the computation of the intrinsic 



shear stress af'^ as a function of the shear rate Kt- Following Chap. 3 in Ref. ITBI or Chap 



16.3 in Ref. |T6|, we express the shear stress in terms of the force per unit area, exerted by 
the polymers 



N 

= ^Xmv^ E E ^"(^' s). (2.8) 

i=l s=l 

Here, Rt(i,s) is the solution of the equation of motion ( |2.2| ) with some initial condition 
Rio(i,s) at time to in the distant past. Moreover, pQ stands for the polymer concentration 
and F"{i,s) := —dH/dRf{i,s) is the force acting on monomer {i,s) at time t. For the 
simple shear flow ( |2.1| ), a linear response relation 

at'^ = I drGt^r^r (2.9) 

with response function Gt will be shown to be valid for arbitrary strengths of the shear rate 
Kt- For a time- independent shear rate k the shear stress is also independent of time and the 
intrinsic shear viscosity r] is then related to the shear stress via 

viG) ■■= fx"'V(«^Po) = Po' / dt Gt. (2.10) 

Jo 



6 



Note that we have made the dependence of t] on the reahzation Q of the crosshnks exphcit. 

Second, in the absence of the shear flow, k = 0, we compute the intermediate incoherent 
scattering function 



Stiq\g):= lira ^ ^ exp{z q ■ (Ri(2, s) - Ro(2, s))} (2.11) 

1 = 1 s = l 

with special emphasis on its long-time behaviour. Due to isotropy the scattering function 
depends only on the modulus q := |q| of the wave vector. 

We expect that observables, like the viscosity and the incoherent scattering function, are 
self averaging in the macroscopic limit and compute the averages 

{ri):=J2P(^)viO) (2.12) 

g 

and, accordingly, {St{q)) over all crosslink realizations. 



C. Formal Solution 

The equation of motion (|2.2|) is linear for the spatially homogeneous shear gradient ( p.l|) 
and hence can be solved exactly. For notational convenience we introduce A^L-dimensional 
vectors such as Rf := {Rf (1,1), Rf{N, L)), whose NL components are the respective 
spatial components of the position vectors Rt(i,s). The force vectors = —{d/a^)TRf 
and the noise vectors , a = x,y, z, . . ., are defined in an analogous manner. Furthermore 
we use the abbreviation K^'^ = Sa,x^/3^yKt for the spatial components of the velocity gradient 
tensor, which is according to ( |2.1| ) spatially homogeneous. 

The expression ( p.8|) for the stress tensor can be rewritten as 

af'^ = ^ lim Tr(rcr'^) 

, N L 

i=l s=l 



where C"'^ := Rf{R^)^ is the matrix of second moments of Rt and the superscript"'" denotes 



transposition. The matrix C"'^ may be calculated from the equation of motion ( |2.2| ) as 
follows: Introducing Utii, s; i', s') := exp{—dt/{(a^) F} (z, s; i', s'), it is readily verified that 

< = Ut-t, ^^S< + /* dt' f/,_. (2. 14) 



is the unique solution of ( p.2| ), with initial condition R^q, provided that T'^f is the solution 



of the differential equation dt T^f = J2y^t'^ '^tf with initial condition T^f = 6a, /s- Since 



n 

is Gaussian white noise with zero mean, we obtain 
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- U2{t~to) J t,to '-'to t,to 

7:7' 

+ ^ fdt'U,^,.,^ E^*7^5''- (2-15) 



Multiplying Eq. (|2.15|) by V and fixing the initial condition at time to ^ ~oo we end up 
with 

hm r Ct'^ = - f dt' V f/2(,_.) ^ r,7 T57. (2. 16) 



In writing down (|2.16|) we have taken advantage of the fact that F is a positive semi-definite 
matrix, implying Iimj_,oo Ut = Eq, where Eq denotes the orthogonal projector onto the null 
space of r. Inserting ( |2.16| ) into ( p.l3| ), we find an expression for the stress tensor 



",/9 
't 



a 

.L^ \dt' ^ ^ j ' 

7 

in terms of the time- dependent linear response function 

2dt 



G. = ^Tr((l-i^o)exp{-^r}). (2.18) 



As a consequence of the simple shear flow (|2.1| ) we have T"^', = 5a,i3 + 5a,xS/3,y j\, dTK{T). It 
follows from definition ( p.lO|) that the static shear viscosity 

is given by the trace of the Moore- Penrose inverseEl of F, i.e. the inverse of F restricted to 
the subspace of non-zero eigenvalues. 

Next we turn to the computation of the incoherent scattering function in the absence of 
the shear flow. In this case we have k = and hence T^f = 6a,f3- Note that Rt — Rq is 
a Gaussian Markov process whose distribution is characterized in the limit tQ ^ —00 hj a 
vanishing mean and the covariance 



to 



lim^(K-i?o) {R?-R'o] 



^6^,^iyTexp!^-^^r}. (2.20) 



To derive ( |2.20D we have used the solution ( ^.14| ). Hence the scattering function is expressed 



in terms of the diagonal matrix elements of the NL x iVL-matrix 



t ^ a"^ ^ — Eq / _ ( dt ^ , . V 
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via 

^ N L 

St{q\g) = — Y.Y.''M-<l' 9t{G\t,s;t,s)} , (2.22) 

i=l s=l 

where, again, the dependence on the reahzation Q has been made exphcit. 



D. Cluster Decomposition 

Each crosshnk reahzation Q defines a random labeUed graph, which can be decomposed 
into its maximal path-wise connected components or clusters 

K 

g=[j^^k. (2.23) 

fc=i 

Here, Afk denotes the k-th cluster with Nk polymers out of a total of K clusters. The 
associated connectivity matrix, also called Kirchhoff matrix or admittance matrix in graph 
theory, see e.g. pp. 54 in Ref. |2^, is of block-diagonal form 

K 

r(^) = 0r(A4). (2.24) 

k=l 

Therefore the viscosity ( |2.19| ) is decomposed into contributions from different clusters ac- 
cording to 

^(^) = E ^ ^(^^) ' (2.25) 

k=l 

where the individual contribution from cluster A/fc is defined by 

.^U.y.^S^^Jl^J^]. ,,26, 



2dNk V r(M 

In the same way the incoherent scattering function may be decomposed into contributions 
from different clusters 

St{q\g) = Y.J^S,{q\Ar,), (2.27) 

k=l 

where the individual contribution from cluster A4 is defined by 



1 

St{qW := ^ E E exp{-g' 9t{^fk\^, s; s)} . (2.28) 
Combining ( |2.27| ), ( |2.28| ) and ( |2.21| ), we express the scattering function as 
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St{q\Q) = X^^expj 

k=l ^ 



qH 



X exp 



1 



d 



1 



r(A4) 



1 — exp 



dt 



Ca2 



r(A4 



(2.29) 



where we have used the representation 



^1 
Eo{G\i, s; i', s') = ^ 



k=l 



(2.30) 



of the orthogonal projector Eq onto the null space of F. Here A/'(i) denotes the cluster 
containing polymer i. Eq. ( |2.30| ) follows from the fact that the null space of F is spanned by 
the vectors which are constant when restricted to any one cluster. Physically this reflects 
that there is no force acting on the centre of mass of any one cluster. Hence, the number of 
zero modes of F is equal to the total number K of clusters. 

For calculating disorder averages it will be advantageous to reorder the sums in (|2.25|) 
and ( p.29| ) by summing first over all clusters consisting of a given number n of polymers and 
subsequently over all "sizes" n. Thus we obtain a decomposition of the average 



K 



Kk=l I n=l 

for an arbitrary real- valued function /. In ( |2.31 ) we have introduced the average 

(/). 

over all clusters of n polymers, each consisting of L monomers. Furthermore 



f2.311 



(2.32) 



I k=l I 



(2.33) 



represents the average number of clusters of size n per polymer. Note that up to now no 
particular crosslink distribution has been specified. 



III. SHEAR VISCOSITY 



A. Relation to networks of random resistors 

The viscosity ( |2.19| ) of randomly crosslinked polymers can be quite generally related to 
the resistance of a random electrical network. Thanks to the cluster decomposition (p.25|) 
it suffices to consider an arbitrary connected cluster A4- To establish this connection we 
identify a bond between two neighbouring monomers on the same polymer as a resistor 
of magnitude /a? and a crosslink between polymers as a resistor of magnitude 1. The 
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resistance measured between any connected pair of vertices (i, s) and (i', s') will be denoted 
by TZ{i,s;i', s'). For an arbitrary connected resistor network A4 the laws of Ohm and 
Kirchhoff leadEil to a relation between TZ and the Moore- Penrose inverse of the connectivity 
matrix 

n{Nk\t, s]i,s) = — — -r- — (z, s; I, s) 

J- {j\lk) 

l-MUk),., , 

+ r(A4) (^'^'^'^^ 



The matrix elements of the projector Eo^Afk) are all equal to {N^L) ^ according to ( p. 301) . 



Hence, when summing ( p.l| ) over all {i, s), {i', s') G Afc, the last term on the right-hand side 
equals — 2 iV^L Tr[£^o(A/fc)(l — -£'o(-Mc))/r(A4)] and thus vanishes. Using this fact and the 
definition ( p.26| ) of ri{Afk), the viscosity of a cluster A4 is seen to be given by a sum over 



resistances 

^(A4) = ^^ E '^mhs■,^',s'). (3.2) 

{i,s),(i',s')eA/-fc 

Together with (p.25|) this constitutes the announced connection between the viscosity of 
a randomly crosslinked polymer melt and a random resistor network. We would like to 
emphasize that this connection relies on the special form ( p.l9|) of the viscosity in the Rouse 
model. 



B. Finite-dimensional percolation: scaling description 

The gelation and the vulcanization transition have been interpreted as a percolation 
transition already by Flory and Stockmayer. Even though the classical theories had 
to be replaced by modern approache^il built on the analogy to critical phenomena, the 
percolation picture survived. It is generally believed that the ensemble of macromolecular 
clusters which are built in the process of gelation has the same statistical connectivity as 
an ensemble of clusters which are constructed in a d-dimensional bond percolation process. 
This identification is supported by experiment, e.g. the measured cluster-size distribution of 
macromolecules is well described by the critical exponents predicted for bond percolation.i'0 

The most intuitive picture of gelation and vulcanization is related to continuum perco- 
lation, which is expected to be in the same universality class as random bond percolation.! 
In continuum percolation a set of points is randomly distributed in (i-dimensional space. 
The points are assumed to be the centres of spheres with randomly distributed radii. Two 
points are said to be connected if their spheres overlap. Correspondingly, if a crosslinking 
agent is added to a dense polymer solution or melt, then we expect that pairs of monomers 
which are close to each other will be crosslinked with high probability, provided the reaction 
is sufficiently fast as compared to the diffusion time of the polymers. 
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The well-known characteristics of percolation clusters will allow us to determine the 
critical behaviour of the shear viscosity as given by (|2.25|) and ( 2.31 ). For that purpose 



we first recall some quantities which are important in the scaling theory of percolation, 
see e.g. Ref. For definiteness we consider bond percolation on the (i- dimensional cubic 
lattice Z*^. Nearest-neighbour bonds are present with probability p := c/{2d) and absent 
with probability 1 — p. The bonds are interpreted as crosslinks, while the lattice points are 
identified with monomers. Here, the probability p is defined such that the average number 
of crosslinks is given by M = cN . A percolating cluster first appears at a critical crosslink 
concentration Ccrit- The fraction Q of lattice points belonging to the infinite network serves 
as an order parameter and vanishes continuously as the percolation transition is approached: 
Q '^■l^^"* (^(, _ CcritY . The correlation function V{r) is proportional to the probability that 
two vertices which belong to the same cluster are a distance r apart. As the transition is 
approached, correlations become increasingly long-ranged as is indicated by a divergence of 
the correlation length ^ 

e^:= J]r2p(r)^^^-(Ccnt-c)-2^ (3.3) 

At the critical point the correlation function shows an algebraic decay 

nr)l=.^J-'^ r-^'-'^^\ (3.4) 

As far as geometric properties of percolation are concerned, the critical behaviour is deter- 
mined by two independent exponents. Here we choose the exponent v of the correlation 
length and the anomalous dimension rj. The other exponents may be expressed in terms of 
those via various scaling relations, such as 

P=^-[d-2 + r,). (3.5) 

It is a well-established fact of the scaling description of percolation that the average number 
of clusters of size n obeys a scaling law 

Tn = ^"^/((Ccrit " c) Tl") (3.6) 

with the scaling function / decaying faster than any polynomial for large arguments and 
approaching a nonzero constant for small arguments. The exponents a and r are related to 
r] and u by 

T = l + - , fT=— -. 3.7 

d + 2-r]' u{d + 2-ri) ^ ^ 

After these preliminaries we turn to the calculation of the averaged viscosity. For the sake 
of simplicity we consider polymers with L = 1, that is, a network of Brownian particles. This 
is relevant for random networks of macromolecules which are generated by polycondensation, 
starting from small units without internal structure. The relation ( |3.2D between the viscosity 
of a cluster and the corresponding resistance then simplifies to 



12 



r/(A4) = ^ E (3.8) 



We remark that the average which was defined in (p.32|) , is an expectation value over 



all lattice animals of size n and iVr„ refers to the average number of percolation clusters of 
size n. 

We make the basic assumption that {i])^ which, according to ( |2.32| ), only tests clusters 
of a given finite size n, does not acquire any irregularity as the critical point is approached 
and that {r])j^\^^^ ,^ ~ with some yet unknown exponent b. We will see below in ( p.22[ ) 
that this holds true with b = 1/2 in the special case where the crosslinks are distributed 
according to mean-field percolation. Furthermore, we remark that this assumption may be 
circumvented in a more sophisticated approach, as will be explained below. 

Using the properties of the scaling function / and (|3.71 ), Eqs. ( 2.25|) and ( p. 31 ) yield for 
the asymptotic behaviour of the viscosity 



oo oo 
71=1 n=l 



Cent - C)~^^''\ (3.9) 

with a critical exponent A{b) := i6z/((i + 2 — 77) — |z/((i — 2 + 7]). Thus, the remaining task 
is the derivation of the exponent b. 

Since the critical behaviour of the viscosity is believed not to be determined solely by 
large-scale geometrical properties, i.e. by the exponents 77 and the appearance of a new 
exponent is expected. From ( p.8| ) and ( p.32| ) the average over clusters of size n is written in 
the form 



k=l i,i'=l 

xn{i;i')), (3.10) 



which yields 



Ad °° I 1 ^ 



^ n=2 \ iA'=l 



' (Ccrit - c)"('-^)"-<^. (3.11) 



To obtain the second equality in (|3.11D we used enumeration invariance of the system, and 



the asymptotic behaviour for c | Ccrit is obtained from Eq. (2.45) in Ref. |T^. The exponent 
(p, which was first introduced in the context of random resistor networks, governJ^^lllil the 
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growth of the resistance TZ{r) between two points on the incipient spanning cluster, which 
are a large spatial distance r apart: TZ{r) ~ r'^/''. Referred to as the crossover resistance 
exponent^'0 is related! to the spectral dimension, according to = vdf{{2/ dg) — 1), 
where df = d — P/u is the Hausdorff-Besicovitch dimension of the incipient spanning cluster. 
Using again the scaling assumptions ( |3.6|) and {ri)^\^^^ . ~ n'', Eq. ( p.llD amounts to 



(3.12) 



n=l 



which, upon comparison with (p.9|) gives A{b + 1) = (2 — 77)1/ + and, after a little algebra, 
the relation 



6 = a0 = (2/4) 
Thus we have derived the scaling relation 

k = (j)-(3 



1. 



(3.13) 



(3.14) 



for the viscosity exponent k. 

As already mentioned, the above scaling assumption is encoded in a more general scaling 



relation. In Refs. and ^ the authors discuss the generating function 

1 



Z{X,iu) :-- 



exp< 



A2 

2 [r + icu 



(3.15) 



of the distribution of the resolvent of F. Amongst other things they show by means of a 
renormalization-group analysis up to second order in e = 6 — d the validity of the scaling 
relation 



Z{X,UJ) (Ccrit - C)^ 

X ^((ccrit-c)-<^A2,(ccrit-c)i/'^^) 



(3.16) 



with an appropriate scaling function z. We use the fact that the system is translationally 
invariant and relate the viscosity to Z via 



2d 



1-En 



lim 



d 



Z{\,u) 



A=0 



d <:-5 

To derive the second equality we have employed the expansion 

1 i „ l-Eo 



(3.17) 



T + iuj u r 



+ 0{iu), cu^O, 



(3.18) 



in ( |3.15| ). From ( |3.17| ) the result ( |3.14j) is readily re-derived. 
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C. Mean-field percolation: exact results 



To go beyond scaling arguments and to compute the viscosity for all crosslink concentra- 
tions, we have to resort to simpler cluster distributions than the one given by bond percola- 
tion in finite dimensional space. We only consider the simplest distribution which is equiva- 
lent to mean-field percolation, as far as static quantities are concerned. All pairs of monomers 
are equally likely to be crosslinked and no correlations between crosslinks are taken into ac- 
count. More precisely, for a function / depending on the realization Q = {if>, s^] i'^,, s'^}^=i of 
the network we define the average by 

(M ^ N ^ \ 

e=l ^ ' ie,i'e=l Se,s'^=l I 

X /({^e,se;^;,<}). (3.19) 

The static shear viscosity can be computed exactly for this simple distribution, making use 
of results in the mathematical literature on random graphs. The calculations are easiest for 
a network of Brownian particles, that is, polymers consisting of just one monomer each (L = 
1). The generalization to more complex molecular units with arbitrary L — not necessarily 
linear — will be shown to be straightforward. It will not change the critical behaviour as 
compared to the case L = 1. 



1. Brownian particles: L = 1 

The statistical properties of clusters generated according to the above distribution, have 
been studied extensively in the theory of random graphs, as developed in Ref. O. Strictly 



speaking, the ensemble of random graphs considered in Ref. 17 comprises random labelled 

' iV ^ 



graphs consisting of vertices, M = cN edges and each of the y Ij j possible graphs 
is equally likely to occur, i.e. double edges or self-loops are suppressed. Obviously this 
ensemble differs from the one considered in the expectation ( p. 191 ), but one expects both to 
have the same properties in the macroscopic limit. This may be understood by comparing 
the number 0{N) of possibilities to realize a self- loop or double edge in Q with the number 
0{N'^) of possibilities to choose a different monomer. 

In the macroscopic limit there are no clusters of macroscopic size for c < Ccrit '■= \i ^^'^ 



all monomers belong to tree clusters without loops, see Thms. 5d,e in Ref. 17. Moreover, 



according to Eq. (2.18) in Ref. |T^, one has 



n-2 

-2c\n 



r„ = — -(2ce-^T (3.20) 
lcn\ 

for the average number of trees of size n per polymer. Since each cluster Mk of size n = 
is expected to be realized as exactly one of the n"~^ labelled trees 7^ of size n, one has 
Sr„ ^J^k,T-n = ^N,„n) and the expectation ( |2.32| ) over clusters of size n is expressed as an 
expectation over trees 
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{f)n = (^f:^^^^r) = (3-21) 

r„ " \ k=i I r„ 

The last step is understood as follows: The assumption that Mk is a tree of size n fixes the 
number of vertices and edges in Mk- Hence, due to the independence of the crosslinks in 
( p.l9| ), the number of graphs, which may be realized within this assumption, depends only 
on the size n, and all trees 7^ can be proven to occur equally likely as a realization of Mk- 

Accordingly, the probability <^A^~^ Yl!k=i ^J^k,T^ ^ cluster A/fc of size n to be realized 

as the tree 7^ is equal to n^~". 

Note further that the resistance between two vertices i and i' in a tree simplifies consid- 
erably, because there is a unique path connecting i and i' , implying that all resistors are in 
series. Therefore the resistance 7^(i; i') is equal to the number of crosslinks connecting i and 
i' . Hence, we refer to Thms. 1,2 in Ref. for computing 

(TZM)^ = {n-2)\Y: 

u=2 ^ >' 



i^i'. (3.22) 
Obviously, one has {n{i,i))^ = 0. Now (^) and imply 



,,—9 ^ ' 



v=2 

Combining (|2.25|) , (|2.31| ) and ( |3.23| ), and expanding the exponential in ( |3.2CI| ), one obtains 



j=2 v=2 ' 



1=0 

To sum up this expression in closed form, we note that 



m 

J](-l)'(« + m- 0-^(7) 
1=0 



a ^, for m = 0, 

-(a(a + l))"\ form = l, (3.25) 
0, otherwise. 
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For m = 0, 1 this is verified by inspection and for m > 2 by differentiating twice Formula 
0.154.5 in Ref. with respect to a. Hence the result of the u and / sums in (^M ) is equal 
to 1/j and we obtain 



Co! 

8dc 



In 



1 



1 - 2c 



2c 



for L 



(3.26) 



for the averaged viscosity for all < c < |. The viscosity exhibits a logarithmic divergence 
at the critical concentration Ccrit = | corresponding to the critical exponent k = 0. This is 
in accordance with the more general result ( p.l4| ) because for d > 6 one hasi /3 = 1 = (/). 
The behaviour ( p.26| ) of the averaged viscosity is displayed in Figure |I] in the Introduction. 



2. Network of polymers: General L 

The result ( |3.26| ) is readily generalized to networks of crosslinked polymers, all having the 
same number L > 1 of monomers and the same geometric structure. The only difference is 
that now there are two different kinds of resistors with magnitudes I'^/a^ and 1, corresponding 
to intra-polymer bonds and crosslinks, respectively. Due to the fact that the monomer labels 
s are distributed independently from the chain labels i in the expectation ( |3.19| ), it follows 



by the same arguments as above that in the macroscopic limit all poljTiiers are connected 
within tree clusters {ie,i'e}- Thus the average ( ^.21| ) over all possible configurations of 
clusters consisting of n polymers generalizes to 



(/), 



[Ut^ E );^E/({^e,.e;C4}). (3.27) 



In ( |3.27|) the average is taken over equally probable trees 7^ = {ie,i'e\e=i ^^^d over indepen- 
dently chosen, equally probable monomer labels. Given i, i' G 7^, the resistance TZ{i, s; i' , s') 
splits up into two additive contributions, see also Fig. The first is the crosslink part 
TZ{i; i'), which is equal to the number of inter-polymer crosslinks, that is resistors of magni- 
tude one, on the unique path on the tree 7^ from i to i'. Let us denote these crosslinks by 
{iek, Se^,]i'e^., s'^f.}'k=i^ ■ The second stems from intra-polymer resistors of magnitude /^/a^. 
Thus, we write 

7^(^, s; z', s') = 7^(^, + ^ E 2(4' ^^.^i) ' (3-28) 

fc=0 

where s'^^ := s, -Se^^. := s' and {P/a'^)Q{a,a') is the resistance measured between 
monomer a and monomer a' on the same polymer. It follows that 

L /'R.{i,i') \ 

E ( E 2(4' ) = (1 + (7^(^, 0) J Ql . (3.29) 

s,s'=l \ k=0 I n 

The latter holds due to the fact that monomer and polymer labels are distributed indepen- 
dently in the average ( p. 27] ), and we have introduced the averaged single-polymer resistance 
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FIG. 2. Contributions to the total resistance in Eq. ( |3.28| ). Crosslinks are depicted as zig-zag 
lines, intra-polymer bonds as straight lines. 



Upon inserting ( p. 2^ ) into (|3l^), thes 

f f \ CP 



(3.30) 



(T,(t' = 1 

these two relations lead to 



(3.31) 



Thus, the averaged viscosity for general L > 1 follows from ( |2.31| ), ( |3.31| ), the averaged 



viscosity ( p.26| ) of the corresponding Brownian-particle case and the relation ^ 



which is valid for all < c < ^. This gives rise to the exact result 



n=l 



riTn 



(a 



8cd \ 0? ^ 
4d 



In 



2c 



2c 



(3.32) 



for the averaged viscosity of a crosslinked poljTiier melt on the sol side. As in the Brownian- 
particle case, it exhibits a logarithmic divergence at the critical crosslink concentration 
Ccrit = \- The result (|3.32|) is universal in the sense that the details of the model only 
affect the pre-factor of the critical divergence. The pre-factor depends on the persistence 
length /, the "extension" of the crosslinks a, as well as on the number L of monomers of a 
polymer and on the geometric structure of a polymer through the averaged single-polymer 
resistance Q^. Let us mention three concrete examples for the latter. First, in the special 
case where the polymers are chain molecules, the averaged single-polymer resistance is given 
=1 1*^ ~ '^'1 ~ (-^^ ~ l)/(3-^)- Accordingly, for long chains the viscosity is 
proportional to L^, as it should be within a Rouse-type model.lii Second, for ring polymers 
Ql = {L"^ — l)/(6-Z^) is half as big as for linear chains. Third, for star polymers one obtains 
Ql = 2{L — Vf' jl? which implies a linear growth of (77) in L for large L. In all three cases 
the scaling of the prefactor with L in (|3.32|) is not altered when passing to the limit a | of 
hard crosslinks. 
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D. Alternative derivation with replicas 



In this section we rederive the result (|3.26|) for the disorder average of the viscosity of 



a network of Brownian particles by means of the replica trick. We consider the mean-field 
distribution of crosslinks, where each crosslink 1 <i < i' < N, occurs independently 

with probability 2c/ N. In many respects this ensemble is equivalent in the macroscopic 
limit to the one with a fixed number of crosslinks, which underlies ( p.26|) and was defined 
by the expectation ( p.l9| ) for L = 1. In particular, Ccrit = | and the expression ( p.20|) for the 



average number r„ of trees of size n per polymer remain valid, as can be inferred from Thms. 
V.7 and V.5 in Ref. |3^. Additionally, Eq. (|3.21| ) continues to hold, because the crosslinks 
are distributed independently and with equal probabilities. Therefore the averaged viscosity 
is given by (|3.26|) for both ensembles. 

In this subsection we focus again on the non-percolating regime c < Ccrit = ^- To leading 
order in 3> 1 the ensemble, where the crosslinks are distributed with probability 2c/ N, 
is characterized by the characteristic function 



(e^(rA)^ = exp J ^ V e^-+^«'>-'-^-'-^»'.- - cN \ (3.33) 




of the connectivity matrix F. Here, A denotes an arbitrary x A^-matrix. Note that, 
contrary to our previous conventions, in this subsection we do not include the limit N ^ oo 
in the average (•). The ensemble ( p.33|) of random matrices is closely related to that in Ref. 



JT], where the density of states of T was investigated with the replica trick; see also Refs. 



and ^ for recent works devoted to this issue. 

The following computation of (77) with replicas has a lot in common with that of the 
density of states in Ref. ^ Therefore we have chosen a less detailed exposition in this 
subsection and only concentrate on the main steps. Compared to Ref. ^ there is one big 
difference, however. For the purpose of (rj) one can solve the resulting saddle-point equation 
exactly by analytic means. 

The basic quantity in the replica calculation is the averaged resolvent or Green function 

Gi^) ■■ = ^Tr ^ ^ 



N^oo2N\ \T + n 



(3.34) 



of F, in terms of which the viscosity (p^.l9D is written as 



hm (r^) = ^ hm f G{n) - lim ^^'^^'^^ 



= ^ lim (cin) - . (3.35) 



d n-^o \ 2Q ^ 

To get the second equality in ( p.35| ), we have used that the dimensionality of the null space 



of F is equal to the total number of clusters and refer to Thms. V.7(ii) and V.12 in Ref. |3G 
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In the following we compute G{Q) with replicas and show that it has a singular contribution 
for f2 — >^ which cancels the second term in (|3.35|) . The regular part of G{Q) determines 
the viscosity. 

Introducing the generating function 



z^{n) := 

X exp 

the resolvent (|3.34|) is expressed as 



N 

n 



\ i,i'=l 



(3.36) 



(3.37) 



Now the standard replica technique is applied to perform the average of the logarithm 
in ( ^.37| ) via InZ^v = lim„^o(-^]v ~ l)/""-- so doing we introduce hatted vectors x := 
) for n-times replicated variables. Then, specifying Ai^ii := — (xj ■ Xj/)/2 in the 
characteristic function (|3.33|) , we get for the disorder average 



((Z^(fi))") 



N 

n 



N 



L (2vr)«/2 



exp 



i=l 



'iE^^ + ^E^-''^-''''^'-^^]- (3-38) 



N 



i,i'=l 



Next we introduce a delta-correlated Gaussian random field ip : M" M with zero mean 
which serves to decouple the double sum over exponentials in the exponent of (|3.38|) . This 
leads to the functional-integral representation 



(3.39) 



where 



F{n,n,ip) 



c + 



dx (tplx))'^ — In 



X exp 



dx 



2c I - 

TT 



i2n 

n/4 



|n/2 



dy^iy) e-(*-^)' 



(3.40) 



The saddle-point method is now applied to evaluate (|3.39| ) in the limit N oo, yield- 

ing {z^inr) 

equation 



N-^oo 



exp{—NF{Q.,n,ip)) where ip is the unique solution of the saddle-point 
S 



Si; 



F(n,n,V') 



0. 



1p = 1p 



(3.41) 



Note that existence and uniqueness of ip are guaranteed, since for crosslink concentrations 
< c < I it is straightforward to prove that F{Q, n, ip) is strictly convex in the argument 
ijj. For the following it will be more advantageous to work with the Gaussian transform 
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2c ( - 

7r 



n/4 



of the solution ip, in terms of which the saddle-point equation ( |3.4lD is written as 

dy e~ 



(3.42) 



= 2c 

and the resolvent ( p.37| ) reads 



(3.43) 



dyy e 



2 -nj/2/2+</'(y) 



G{n) = lim — 
Note that ( 3.43|) implies the normalization condition 



(3.44) 



dx (f){x) = 2c (27r)"/2 . 



(3.45) 



In order to compute the viscosity ( |3.35|) , the equations ( p.43|) and (|3.44| ) have to be solved 
for small Q after a continuation to n = 0. We restrict ourselves to solutions of ( ^.43p , 
which preserve rotational invariance in replica space and make the replica-symmetric ansatz 



Q \x\ 



(3.46) 



For later purposes the dependence on parameters has been made explicit. Introducing n- 
dimensional spherical coordinates, Eq. ( ^.44| ) transforms into 



G{n) = lim 



dp p"+^e 



dp p"-'e 



n.-l -p2/2+vS(p) 



(3.47) 



To proceed with Eq. (|3.43| ) we employ an integral representation of the exponential of the 
ra-dimensional Laplacian. More precisely, for any rotationally invariant function f{\y\) the 
identity 



n f d^ n-l d 



/(P) 



(3.48) 



p=\x\ 



is valid for all n G N and its application to ( |3.43D yields 
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(f_ n-l d\ 1 2/2+^n(^) 

(3.49) 



exp< 


\VL f n-l dV] 
[2 p dp)\ 


^ g-p2/2+^5?(p) 


'"J 


poo 

i dr] 77"-^ e"'''/' 







Here s„ := n7r"'/^/r(n/2 + 1) denotes the surface of the unit sphere in MJ^. The normahzation 
( p.45| ) translates to 

dp p"" Vn (P) = 2c (27r^])"/2 . (3.50) 

Observing, that hm„^o /o°° '^P P"~^/(p)) = /(O) is vahd for any sufficiently fast decaying 
function /, we get from ( p.49| ) and ( p. 47] ) 

.?(p,^2ee-exp{§(|,-iA)}eV/-i>. (,51) 



G{n) = ^e-^'' I dppe-p'/^+^oip) ^ (3.52) 



and 

2n 

respectively. Here we made use of 

ip^{0) = 2c (3.53) 

which is obtained from ( p.50|) . The remaining task is to compute ^q{p) up to first order in 
Q, which, in turn, determines the 0(l)-contribution of G{Q). 
To this end we make the expansion 

^^{p)=V>oip)+9ip)^ + 0{n'), (3.54) 
insert it into ( p.51j ) and equate the zeroth- and first-order terms 

»'» = -W'(-2ce-2'-''V2) , (3,55) 



(3.56) 



Here, W denotes the principal branch of the Lambert Vr-function,Ej which is defined as the 
real solution of the transcendental equation 

W{x) exp{W{x)) = X, x>-l/e. (3.57) 

Moreover, from inserting (p.54|) into ( p.52| ), we get the desired expansion 

1 r°° 

G{n) = -j^ dpp4{p)[Q-'+g{p)] +0{n). (3.58) 
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In order to calculate the integrals in (|3.58|) , we employ the relation 

dp 1 - foip) ' 

which follows from ( |3.55|) and ( p.57|) . It helps to rewrite both the first term 



P^oip) 



d 
dp 



4ip) { 2 Voip) - 1 



(3.59) 



(3.60) 



and, in combination with ( p.56| ), the second term 
P<^o(p)£/(p) 



1 d 
'2dp 



d 

dp 



^°(p) - In (1 - ^°(p)) - ^°(p) 



(3.61) 



of the integrand in (|3.58|) . Observing in addition the normalization (|3.53|) and {dipQ/dp){0) 
0, which also follows from ( p.59| ), the evaluation of the integral in ( p.58| ) leads to 



Gin) = ^ - ^ [ln(l - 2c) + 2c] + 0{n) 



(3.62) 



Hence, the singular contributions for Q 
expression 

lim (ri) = - — - 



In 



cancel in (|3.35|) , and we end up with the 
1 



1 - 2c 



2c 



(3.63) 



for the averaged viscosity, which coincides with the exact result ( p.26|) derived in the previous 
section. 

We conclude that the replica trick and the assumption of replica symmetry provide an 
exact method to perform the disorder average in the random matrix ensemble ( |3.33|) below 
the critical concentration. This was demonstrated for the viscosity, but we conjecture that 
it holds more generally, e.g. for the density of states. 



IV. THE INTERMEDIATE INCOHERENT SCATTERING FUNCTION 

Dynamic density fluctuations in polymer melts are encoded in the intermediate incoher- 
ent scattering function ( p. 11 ). The leading long-time behaviour of its average {St{q)) over 



all crosslink realizations was computed in Ref. ^ for the ensemble ( |3.19| ) of uncorrelated 
crosslinks, for which Ccrit = |- Among others, it was found that 

(1) the incoherent scattering function decays for all crosslink concentrations < c < Ccrit 
like a stretched exponential 

{Stiq)) * t-'/' exp{-it/t,r} with a = ^, (4.1) 
where tg ~ is a diffusive time scale. 
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(2) as the gelation transition is approached the diffusive time scale diverges 

q\ '^^^'^ (c„it - c)-^ with /i = 2 , (4.2) 
and the effective diffusion constant vanishes linearly 

(Ccnt-c)-^ (4.3) 



eff 



where y = 1 and Dq := (CL) ^ is the bare diffusion constant, 
(3) the incoherent scattering function decays algebraically at the critical point 

{St{q)) * with x=^. (4.4) 



The analysis in Ref. |^ relies on the cluster decomposition ( p. 271) for St{q) in the form (|2.31|) 



{S,{q)) = J2nTn{S,{q))^. (4.5) 



n=0 



In addition, the plausible assumption has been used that the long-time asymptotics of {St{q)) 
is obtained by inserting the long-time asymptotics of {St{q))^ into the right-hand side of 

(PD- 

The purpose of this section is twofold. First, for the case of uncorrelated crosslinks as 
in Ref. ^, we check whether the corrections to the leading long-time behaviour of {St{q))^ 
do not influence the leading long-time behaviour of {St{q)), as was assumed in Ref. ^. This 
is done in Subsection [lV.lj| by asymptotic evaluations of an upper and a lower bound on 
{St{q)), which are constructed in Subsection |IV.A| . Second, in Subsection |IV.Q we determine 



the exponents a, ^, y and x from the known critical exponents of random bond percolation 
by scaling relations. This generalizes the results in Ref. |] to a non-mean-field distribution 
of crosslinks. 



A. Upper and lower bound on {St{q)) 
Since the connectivity matrix F is positive semi-definite, the second exponential in (|2.29D 



is bounded above by 1, hence 

(4.6) 



Stiq) < X^^expl- 
k=i ^ 



For the computation of the crosslink average it is advantageous to reorder the sum according 
to ( p3TD 
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{St{q)) <Ut{q) :=f;nx„exp|-:^| . (4.7) 

n=l n ) 

A lower bound on St{q) is obtained by neglecting the double-exponential contribution in 



Next we apply the Jensen inequality to the right-hand side of ( |4.8| ) by replacing the nor- 
malized i- and s-sums over the exponentials by the exponential of the normalized sums. 
Together with (|2.31| ) this yields 



{St{q)) > ^?2r„exp 



n=l 



DoqH 



n 



The lower bound is completed by yet another application of the Jensen inequality, 
(exp(«))^ > exp((«)^), and the identification (|2.26|) of the resulting exponent with the 
viscosity 

(5t(g))> A(g):=fjnr„exp|-Dog'(^^ + 2(r7)„^| . (4.10) 

In summary, we have found the chain of inequalities 

Ct{q)<{Stiq))<Ut{q). (4.11) 



B. Uncorrelated crosslinks 

It was argued in Ref. || that {St{q))^ exp{—Doq'^t/n}. This relation follows from 
( ^.29| ). Assuming that the corrections to this leading behaviour do not play a role in the 
cluster decomposition ([4.5[), the Kohlrausch law 



{Stiq)) ^nTnexp{-Doqh/n} 

n=0 



(4.12) 



< c < |, was found! which yields the critical exponents in ( [4.1|) - ( |4.4| ). Here we introduced 
the function h{c) := 2c — 1 — In (2c) and observed h{c) ~ (ccrit — c)^ for c jcrit- In order to 
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derive the second line in ([4.12|) , the exphcit form ( |3.20|) of Tn in the case of uncorrelated 
crosshnks was used, n\ was replaced by its Stirling asymptotics, and the sum over n by an 
integral, whose evaluation yielded the result. 

In this subsection we want to confirm that the corrections to the asymptotic behaviour of 
{St{<l))n ^'^ ^'^^ influence the Kohlrausch law (|4.12| ). To do so we employ the bounds ( [4.111 ). 
We also perform a more careful evaluation of the sums over the cluster sizes n, than the one 
which was used in Ref. ^ to derive the second line in ( |4.12| ). Details of this calculation are 
deferred to the Appendix. 

We start with the upper bound ( [4.7D which reads in the notation (ST) of the Appendix 
as Ut{q) = F{0, D^q^t). Therefore we conclude from (|A19| ) that 



hn. ^ 

et^oo {8cWoqH)-y^ exp{-2[h{c)DoqH]y^} 



1. 



(4.13) 



This result also confirms the previously obtained asymptotic equality in the second line of 

As to the lower bound ( [4.1UD we first recall the expression ( |3.31| ) for (rj)^ and the in- 
equality {r])^ \^^^ < by/n, which follows from (|3.23|) and (|3.22| ), with b > being a finite 
constant. This yields, again in the notation (|A1|) of the Appendix, 



where Bi 
implies 



Ct{q)>e-''^^'F{q'B2,Doqh), (4.14) 
DoCPLQL/{2d) and B2 := DobCa^L{l + QLp/a^)/{2d). Consequently, (|A20|) 



lim 

q'^t — >oo 
q^t->0 



A(g) 



^cWoqH)-y^exp{-2[h{c)DoqHY/^} 
Taken together, ( [4.13|) and (|4.15|) provide us with the equality 

ito 

(8c2B„g2i)-i/2exp{-2|ft(c)£>„i,2i]i/2} 



> 1. 



(4.15) 



1. 



(4.16) 



In other words, the bounds (|4.11| ) only confirm the previously obtained Kohlrausch law 
( [4.12| ), if one also considers the limit of small momentum transfer — as may be a reasonable 
approximation for an experimental low-angle-scattering situation. Nevertheless, we believe 
that the stretched exponential ( |4.12| ) reflects the true long-time asymptotics of {St{q)) for 
finite q, too, and that the lower bound ([4.10|) is not sharp enough to reproduce this. 
On the other hand the bounds ( [1.1 1| ) are sharp enough to establish the expression 



00 

-1 ^ ^ 



n=0 



Do{l-2c) 



(4.17) 



for the effective diffusion constant D^s- The first equality was already stated in ( [1.3[) . The 



second follows for < c < | from differentiating the relation 1 = Yl'^=i ''^^n with respect 
to the crosslink concentration c. Obviously, the result is in accordance with y = 1 in ( [4.3[) . 
The linear dependence of D^s on c is displayed in Figure |1] in the Introduction. 
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C. Scaling description 



In this subsection we extend the results of Ref. | to the crosslink statistics generated by 
d-dimensional random bond percolation, whose basic properties were recalled in Subsection 
TTTBI. 



For the case of uncorrelated crosslinks we have seen in Ref. |g — and confirmed in the 
preceeding subsection by alternative methods — that the critical exponents ( [4. ID - (O) of 
the Kohlrausch law ( ^.121 ) do not depend on the number L of monomers per polymer. 
Therefore we restrict ourselves to the case L = 1 of Brownian particles in the following. The 
results of the last subsection also motivate us to replace {St{q))n by its long-time asymptotics 
in (O), that is, 



{SM)) 5^nr„exp{-Z}og't/ri} • (4.18) 

n=0 

For t ^ oo this series is expected to be dominated by terms with large n. Therefore we insert 



the scaling law (|3.6|) for r„ in ( 4.18 ) and use the fact that the associated scaling function / 
decays exponentially for large n. This yields up to a constant 



oo 



{St{q)) 5^ni-^exp{-(c,Ht - cf/'^n - D,qH/n] . (4.19) 

n=0 

In order to proceed, we replace the sum over n in ( [4.19| ) by the integral 

POD 

/ dn r/-^ exp{-(ccrit - cf'^n - D^qh/n} , (4.20) 
Jo 

which should not affect the long-time asymptotics either. Using Formula 3.471.9 in Ref. |2^, 
we then express the integral in ( [4.20|) in terms of a Bessel function. Its asymptotic expansion 
according to Formula 8.451.6 in Ref. ^finally gives, up to a constant, 

^ {D,qHY-^-^)l^ 

X exp{-2[Dog'(ccnt - cf'^tYl'} (4.21) 
for c < Ccrit- For c = Ccrit we conclude 

directly from ( [4.201 ). Hence, we find a = 1/2 for the critical exponent ( [4.1D , the same value 
as was found for uncorrelated crosslinks. In contrast, the exponents defined by the relations 
( |4.2|) - (|4.4|) turn out to be different. More precisely, we get the scaling relations 

/i = l/o- ^ 2.22, (4.23) 
1/ = (3 -r)/(T ^ 1.82, (4.24) 
a; = r-2^0.18. (4.25) 

The approximate numerical values are those for random bond percolation in three 
dimensions.0 A discussion of these results in comparison to the experimental data will be 
given in the next section. 



27 



V. CONCLUDING REMARKS 



Starting from a Rouse model of a crosslinked polymer melt, we discuss the critical dy- 
namics of the gelation transition with particular emphasis on the static shear viscosity and 
the long-time behaviour of the incoherent scattering function. Our main result is an exact 
relation for the critical exponent of the viscosity, 

k = (j)-P, (5.1) 

in terms of the crossover and thermal exponents of percolation theory. Two crosslink distri- 
butions have been analyzed in detail: one corresponding to a mean-field percolation model 
and one corresponding to finite-dimensional percolation. For the first distribution we were 



able to derive the exact expression (|3.32|) for the static shear viscosity. This result is valid 
for all crosslink concentrations < c < Ccrit = | and exhibits a logarithmic divergence at the 
critical crosslink concentration Ccrit, that is, k = 0. The critical exponent k = also follows 



from the scaling relation (|5.1| ) when inserting the mean- field values = 1, see e.g. Ref. |23, 
and (3 = 1. For the second crosslink distribution the scaling relation ( |5.1| ) can either be eval- 
uated in terms of the e = 6 — d expansionil for cj) and P, yielding k = e/6 + lle^/1764 + 0(e^). 
Alternatively, one may insert the numerical value of (p known from high-precision simulations 
of random resistor networks in c? = 3 dimensions,!! which, together with the corresponding 
value for (3 gives k\d=3 ~ 0.71. 

The experimental values for k are systematically larger. Early experiments by Adam 
et a/.i were performed near the gelation threshold of polycondensation reactions. For the 
samples with low molecular weight, data for the viscosity were fitted to a power-law diver- 
gence with exponent k = 0.8 ± 0.1. On the other hand, the viscosity data of samples with 
high molecular weight could not be fitted to a power-law divergence. Considerably larger 
exponent values were observed in more recent experiments on silica gels and epoxy raisins. 
Colby et alM report values 1.4 < k < 1.7 and Martin et al. obtained k = 1.4 it 0.2 from 
viscoelastic measurementsU and k = 1.5 ±0.2 rather indirectly from measurements of the in- 
coherent scattering function.i The origin of the scatter of the experimental data is not clear. 
A splitting of the static universality class into different dynamic ones has been suggested.i 
Another possible explanation is the size of the critical region, which is expected to depend 
on chain length.il The observed range of exponent values could then be interpreted as a 
crossover from mean-field to critical behaviour. 

What are the shortcomings of our theory? We have written the average macroscopic 
viscosity as a weighted sum of contributions from individual clusters. Such a decomposition 
holds exactly within the Rouse model of a crosslinked melt, but is expected to be valid 
more generally — as long as there are no interactions between different clusters. Retaining 
the cluster decomposition, the only unknown is the contribution of an individual cluster, 
because the distribution of clusters should indeed be given by percolation theory and hence 
is known to high precision. As far as single clusters are concerned, we expect that relaxation 
times are longest for Rouse dynamics. If, for example, we were to use Zimm dynamics 
together with a pre-averaging approximation, we would find shorter relaxation times as 
compared to the Rouse model and hence even smaller values of k. This suggests that 
the discrepancy between theory and experiment cannot be cured, if we retain the cluster 
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decompostion. Hydro dynamic interactions, excluded-volume interactions and entanglement 
effects, all invalidate the cluster decomposition. One expects entanglement effects to play 
a vital role in stress relaxation.El However, the static viscosity measures stress relaxation 
only on the longest time scales in the sol phase. It has been arguedHS that in the regime 
close to the transition, entanglement effects are less important because of two reasons: first, 
there are almost no permanent entanglements in the sense of interlocking loops. Second, the 
time scale of a temporary entanglement of two clusters is determined by the smaller cluster, 
whereas the dynamics on the longest time scales is determined by the larger cluster. From 
calculations of static quantities we know that the excluded-volume interaction is essential 
for the stability of the network in the gel phased In the sol phase, on the other hand, the 
role of the excluded-volume interaction is less transparent. On one hand, the e = 6 — d 
expansion for a gelation model with excluded-volume interaction yields the same critical 
exponents as obtained from the e-expansion of percolation theory.! On the other hand, the 
excluded-volume interaction is known to be relevant^ for static and dynamic properties of 
polymer melts and solutions below d = 4. 

The decay of density fluctuations in silica gels has been investigated by quasi elastic light 
scattering.!! In the sol phase, a stretched exponential of the autocorrelation was observed, 
St{(l) ~ 6xp{ — (t/tq)"} with a = 0.65 ±0.05. The experimentally determined timescale tg is 
diffusive and diverges as the gelation transition is approached tg ~ (ccrit — c)~^'^. The critical 
behaviour of the diffusion constant was determined as well and in particular the exponent 
value y = 1.9 ± 0.1 was found. At the critical point, the scattering experiments reveal an 
algebraic decay in time of St{q) with an exponent x = 0.135 ± 0.015. All these findings 
are in qualitative agreement with our theory. In fact our expression ( |2.29| ) for St{q) has 



been suggested on phenomenological grounds as a starting point for the discussion of the 
critical dynamics at the sol-gel transition.0 The exponent values of mean-field theory, see 
Eqs. ( [4.1|) - (f4.4| ), deviate from the experimental ones, as one would expect. If we use the 
scaling description of finite-dimensional percolation, then the exponent a of the stretched 
exponential is unchanged. Its value, a = 1/2, is characteristic of Rouse dynamics and 
independent of cluster statistics. The corrections to the other exponents go in the right 
direction, see Eqs. ( |4.4| ). However, discrepancies beyond the experimental error bars remain 
and are probably due to our neglect of excluded-volume and hydrodynamic interactions. 

Our approach can be extended to study stress relaxation at finite frequencies. Within 
the Rouse model, the dynamics is completely determined by the eigenvalue spectrum of 
the random connectivity matrix F. For the mean- field distribution of crosslinks the eigen- 
value spectrum can be computed with the replica trick0 so that stress relaxation at finite 
frequencies becomes accessible to analytical theory. Work along these lines is in progress. 
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APPENDIX: EVALUATION OF THE BOUNDS 

Here we present some technical details which are needed in the asymptotic evaluation 
for t —>■ oo of the bounds ( [4.11|) on the averaged intermediate incoherent scattering function 



St{q). 

The quantity of interest is the series 



F{B, T) ■= J2 e-^/"e-^v^ , (Al) 



n=l 



where 5 > 0, T > and r„ = n"-2e-"^(^V(2cn! e'^) follows from with h{c 



2c — 1 — ln(2c) . Now, pick a natural number P, split the series in two parts given by the first 
P — 1 terms and the terms with n > P, respectively, and apply the Stirling approximation, 
Eq. (6.1.58) in Ref. |g, 

n\ = (27r)i/2 ^n+i/2 exp{-n + 9{n) /{I2n)} , (A2) 

< 6{n) < 1, to the terms with n > P. Thus, we infer the existence of a constant 

e-i/(i2P) < < 1 (A3) 



such that 



where 



FiB, T) = RpiB, T) + ^^^^ F{B, T) , (A4) 



Rp{B,T) ■.= Y,{ 

n=l ^ 



P-1 

" 2c (27r) 1/2^3/2'' 

-T/n ^-B^ 



and 



xe-'/"e-^^" (A5) 



F{B, T) ■= n~^'^ exp{-n/i(c) - T/n - B^} . (A6) 



n=l 



Note that 



^-1 -nh{c) 



n=l 



decays exponentially for T oo, uniformly in i? > 0. 

In order to proceed with F{B,T) we use the Fourier representation 



-xV(2a) 

e-^/^ = dx e'"" - — a>0, (A8) 



(27ra)V2 ' 
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and arrive at 



F{B,T) 



dx 



{An Doty 



poo ^ 
■^M^) n=l 



(A9) 



where A{x) := h{c) +xV(4T). According to Eq. (11.1.1) in Ref. || the series in ( |A9| ) admits 
an integral representation such that 



F{B,T) 



B 



dx e* 



d\ 



A{x) 



di 







(AlO) 



After a partial integration with respect to x and the subsequent changes-of- variables z :- 
a;(4T)-i/2 and C := (20"^/^ Eq. ( [ATq ) reads 



F{B,T) 



<V2 



dc 



(7rT)V2 - (2^)1/: 



dz 



z e 



^z^+hicB/Q _ 1 



(All) 



with h{c,r) := h{c) + The 2;-integration can be done with the help of the residue 

theorem. To this end we close the integration contour with a semicircle in the complex upper 
half-plane and remark that the integrand has simple poles Zn{B/() := z'^{B/() + iz'^{B/(), 
n E Tj, in this half-plane whose real and imaginary parts are given by 



sgn(nj 



\l{h{c,r)f + (27rn)2 - h{c, r) 
^J{h{c,r)f +{27iny + h{c, r) 



1/2 



(A12) 



Hence, we find 



FiB,T) 



1/2 



dC 



-2[Th(c,B/C)]i/2 



(27r)V2 
+2Fi(5/C,T) 



where 



Fi(r,T) := ^e"'^^'')^ cos[2< (r)VT] . 



(A13) 



(A14) 



n=l 



Using z'^{r) > z"{0), we get the r-independent upper bound 



|^i(r,T)| <e-2^"(°)^5^e 

n=l 

< jg-24'(0)v^_ 



-2K(0)-4'(0)]VT 



(A15) 
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Here, / is a uniform constant for all T > Tq > 0, because z"{0) > z'KO) for n > 1 and 
<(0) ~ forn ^ oo. Thus we conclude from (^) and (|M5D that 



y FiO,T) 

(8c2r)-V2exp{-2[/i(c)T]V2} 



Up 



(A16) 



holds for all P G N. 

For B > the inequality ^/h{c, r) < \/h{c) + 2~^/2|r| leads to the lower bound 



(27r)i/2 ^ 

> g-2[T/i(c)]i/2 



(27r)V2 



(A17) 



This bound, together with the same arguments which led to ( |A16|) , now yield 

F{B,T) 



lim 



> Up 



(A18) 



T-^^^ (8c2T)-V2exp{-2[/i(c)T]V2} 
for all P G N. Letting P ^ oo and observing (|A3|), Eq. ( A16| ), respectively ( A18 ), implies 

P(o,r) 



(8c2T)-V2exp{-2[/i(c)T]V2} 



(A19) 



respectively 



lim 



F{B,T) 



T^^^^ (8c2T)-V2exp{-2[/.(c)T]V2} 



> 1 



(A20) 
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